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Abstract. Consider a fluid composed of two species of particles, where the 
interparticle pair potentials u\x = U22 ^ Ui%. On confining an equal number of particles 
from each species in a cavity, one finds that the average one body density profiles 
of each species are constrained to be exactly the same due to the symmetry, when 
both external cavity potentials are the same. For a binary fluid of Brownian particles 
interacting via repulsive Gaussian pair potentials that exhibits phase separation, we 
study the dynamics of the fluid one body density profiles on breaking the symmetry of 
the external potentials, using the dynamical density functional theory of Marconi and 
Tarazona [J. Chem. Phys., 110, 8032 (1999)]. On breaking the symmetry we see that 
the fluid one body density profiles can then show the phase separation that is present. 

In equilibrium density functional theory (DFT) pQ, the key quantity is the 
Helmholtz free energy functional. Given this functional, one can calculate the (ensemble) 
average one body density profile(s) for a fluid of particles subject to any given external 
potential. The form of the Helmholtz free energy functional depends, of course, upon 
the particular interaction potentials between the particles in the fluid. However, this 
functional is in general unknown, except in a few cases for 1 -dimensional fluids pQ. 
None the less, in the last 30 years, a number of rather accurate approximate Helmholtz 
free energy functionals have been developed for a wide class of fluids - see Ref. 
and references therein. When attempting to construct a theory for the dynamics of 
the (ensemble) average density profile of a fluid confined in a time dependent external 
potential, it is very appealing, given this large body of work, to base a theory for the 
dynamics upon the equilibrium Helmholtz free energy functional. In other words, if one 
is wanting to study the dynamics of a fluid, is it possible to just apply the relevant 
equilibrium Helmholtz free energy functional in the dynamical theory? 

Recently, Marconi and Tarazona derived a dynamical density functional theory 
(DDFT) [21 E] for fluids of Brownian particles. Their theory is based upon the 
equilibrium fluid Helmholtz free energy functional in a rather simple way and therefore 
has attracted much interest. The theory is not exact; it involves approximating the two 
body correlations in the non-equilibrium fluid by those of an equilibrium fluid with the 
same one body density profile El EJ However, for all the cases where comparison 
with Brownian dynamics (BD) simulations has been made, the theory has proved to 
be remarkably accurate j2j El El 13 Ej • It appears that as long as one has an accurate 
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approximation for the equilibrium Helmholtz free energy functional, then the DDFT 
gives a good account of the dynamics for Brownian fluids. 

The present paper builds upon an earlier study [Ij, in which we considered two 
different asymmetric binary mixtures of particles interacting via repulsive Gaussian 
pair potentials [Oj: U{j(r) = e^- exp(— r 2 /i?^), where i,j = 1,2 label the two different 
species of particles, e^ > is a parameter that determines the strength and Rij denotes 
the range of the interaction potential. For different choices of these parameters one can 
obtain either a mixture exhibiting bulk fluid-fluid (macro)phase separation [TU I ITT | I7j. 
or alternatively, a mixture exhibiting microphase separation ^21 EZj • We found that the 
DDFT used in conjunction with an extremely simple approximate (RPA) Helmholtz free 
energy functional was able to account accurately for the dynamics of phase separation 
induced by reducing the radius of the confining spherical cavity potential 0. In the 
present paper we consider similar situations but for symmetric binary Gaussian core 
model (GCM) fluids, i.e. where Rn = R 22 and en = e 22 , but where R 12 ^ Rn or 
£i2 7^ £n- The presence of this symmetry can have a dramatic effect on the confined 
fluid one body density profiles in cases where there are equal numbers of particles for 
each species. As long as the confining cavity potentials are the same for both species 
of particles, then the ensemble (or time) average one body density profiles of the two 
species must be exactly the same. Symmetry dictates this must always be the case, 
even when the fluid exhibits phase separation. Here we show that when one breaks the 
symmetry of the confining cavity potential, such that the cavity walls favour one of the 
species, then the fluid one body density profiles can change significantly and show the 
influence of the phase separation - see for example the results in Figs. 

We now describe further the model fluid and the DDFT: The GCM fluid pair 
potentials have no hard core and the centres of the particles can overlap completely. 
Such potentials arise when considers the effective potential between the centres of mass 
of polymers in a good solvent, see Ref. jHj and references therein. In this case ~ 2ksT 
and R^ ~ R g , the polymer radius of gyration. We consider two symmetric binary 
mixtures of GCM particles. In both cases en = e 22 = ei 2 = 2/c#T and Rn = R 22 . The 
first mixture, which we denote "fluid A", is that with i? 12 = 0.6i?n- Fluid A exhibits 
microphase separation, and is similar to the fluid described in Ref. ^2]. This mixture 
exhibits 1-2 ordering, i.e. where a particle of one species preferentially has particles of 
the opposite species as nearest neighbours |T2"] . The second mixture, which we denote 
"fluid B", is that with Ri 2 = l.lRu. Fluid B exhibits macrophase or bulk phase 
separation. This is driven by the fact that R\ 2 > Rn = R 22 ^JJ . The bulk critical point 
for this mixture is at a total density pR\i = 0.54 and, by symmetry, at concentration 
X\ = x 2 = 0.5. 

When the GCM fluid density becomes sufficiently high, the following mean-field 
Helmholtz free energy functional is rather accurate [H| ITIH ITT]: 

F[{Pi(r)}\ = F id [{Pi(r)}]+lY. J dr J dr' A (r) /0 ,(r')^(ir-r'|)+^ J drV^v) Pi (v) , (1) 
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Figure 1. Fluid A density profiles Pi(r,t) (solid line: DDFT results for species 1, 
dashed line: species 2; symbols are BD results, (+) for species 1, (x) for species 2) 
for Ni = N2 = 100 particles, which is initially (t < 0) at equilibrium confined in 
(asymmetric) external potentials with E\ = lOfcsT, E2 = 20/cbT and 1Z = 3i?n. At 
t = 0, E2 is suddenly reduced to E2 = lOfc^T, giving symmetrical confinement. The 
profiles are plotted for various t* = kBTTR^t. This model fluid exhibits microphase- 
separation: the initial configuration has an 'onion' structure, but on enforcing the 
symmetry in the external potentials, any signature of this structure in the one body 
density profiles disappears as time progresses. 



where {pi(r)} are the fluid one body density profiles, Vi(r) is the external potential for 
particles of species i and F id is the ideal gas contribution to the free energy [T|. 

For a fluid of N = Ni + N 2 Brownian (colloidal) particles, one can approximate 
the equations of motion using the following stochastic Langevin equations of motion: 
r^ 1 dr„(t)/dt = —V n U(r N ,t) + G n (t), where r n is the position of the n th particle, Tr 1 
is a friction constant for particles of species i (we assume Ti = T 2 = T), U(r N ,t) is the 
potential energy and G n (t) is a stochastic white noise term - see Refs. [21111 EI for more 
details. From such equations of motion, averaging over all realisations of the stochastic 
noise, one can obtain the following equations for the time evolution of the ensemble 
average one body density profile [21 EH El E] : 
dpi{r,t) 
dt 

where the Helmholtz free energy functional F[{pi(r,t)}] is given by the equilibrium 
functional with the equilibrium density profiles {p«(r)} replaced by the set of non- 
equilibrium profiles {pi(r,t)}. We solve the DDFT for the binary GCM confined in 
spherically symmetric external potentials of the form OEj: V{(r) = (r/TZ) 10 , where 
r is the distance from the origin, Ei is an energy scale and the length-scale 1Z (cavity 
radius) is the same for both species of particles. 

In Fig. [T] we display the results for N\ = N 2 = 100 particles of each species of a 
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Figure 2. Fluid A density profiles Pi(r,t) (solid line: species 1, dashed line: species 
2) for Ni = N2 = 6000 particles. Initially (t < 0) the fluid is at equilibrium confined 
in symmetric external potentials with E2 = E\ = lOfc^T and 1Z = 8i?n. At t = 0, E2 
is suddenly increased to E2 = 20ksT, breaking the symmetry. The profiles are plotted 
for various t* — ksTTR^t. Initially there was no signature of microphase separation 
in the density profiles, but after breaking the symmetry, the fluid profiles develop an 
'onion' structure. 



mixture of fluid A. Initially (t < 0) the external potential parameters E 2 7^ E\ and the 
fluid one body density profiles exhibit 1-2 ordering (microphase separation) in the cavity. 
Then at t = the external potentials are made symmetric: E\ = E 2 . At equilibrium, for 
such symmetric external potentials, the one body density profiles for each species must 
be the same since Aq = N 2 . We see in Fig. Uthat the density profiles quickly evolve 
to give the same final profiles for each species. We also display the results from BD 
simulations, which show extremely good agreement with the results from the DDFT, 
giving confidence in the DDFT approximation. In Fig. |2] we display DDFT results for 
a much larger number of particles. The external potentials for t < are symmetric and 
therefore the ensemble average density profiles at t — are the same for both species. At 
t = the symmetry in the external potentials is broken, such that particles of species 
1 are favoured by the walls of the cavity. The one body density profiles then evolve 
towards equilibrium density profiles that display the microphase ordering in the fluid 
|12j . Note that this does not mean that there was no microphase separation for t < 
(the average total density and composition has not changed); if one were to look at a 
'snapshot' of the particle configurations at a time t < one would indeed see microphase 
ordering. We illustrate this below for fluid B. 

In Fig. El we display the results for Aq = N 2 = 100 particles for fluid B. Initially 
(t < 0) the external potential parameters E 2 = E\ and the fluid one body density 
profiles are the same for each species. At t — the symmetry in the external potentials 
is broken: E\ 7^ E 2 , such that the cavity walls now favour species 1. The equilibrium 
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Figure 3. As in Fig. □ but for fluid B with Ni = N 2 = 100. The initial density 
profiles are those for symmetric external potentials: E\ = E2 = lOfc^T. At t = 0, 
E2 — » 20fcsT. The subsequent profiles then show the signature of phase separation. 



density profiles for the asymmetric external potentials then show that the fluid does 
indeed exhibit the signature of bulk phase separation. We also display the results from 
BD simulations, which show extremely good agreement with the results from the DDFT. 
In Fig. |3] we display the DDFT results for fluid B with a larger number of particles. Once 
again the external potentials for t < are symmetric and therefore the ensemble average 
density profiles at t = are the same for both species. At t = the symmetry is broken 
and the one body density profiles then evolve towards ones that show the fluid-fluid 
phase separation. Interestingly, this evolution to the final equilibrium density profiles 
proceeds via an intermediate state where the fluid is rich in species 1 at both the centre 
and close to the wall of the cavity. Then, as time proceeds, species 1 particles at the 
centre of the cavity diffuse to the outside, since the intermediate state has two fluid-fluid 
interfaces which has a higher free energy cost than the final equilibrium configuration 
with only one fluid-fluid interface. Despite the fact that the t = profiles show no sign 
of phase separation, in the final frame of Fig. 0] we see that in a particular configuration 
of the particles that there is phase separation in the cavity - i.e. the t = profiles 
in Fig. |3] are the result of ensemble averaging over strongly asymmetric configurations. 
This implies that in this larger cavity there are strong fluctuation effects over the whole 
cavity, which are neglected in the present mean field DFT treatment [H|. We therefore 
expect the DDFT results of Fig. 0] to be less reliable than those of Fig. 

We conclude by emphasising that although the initial (t = 0) one body density 
profiles in Figs. |2] El show no obvious sign of phase separation, this does not mean that 
there will be no signature of phase separation in a particular equilibrium configuration 
of the particles subject to symmetric external potentials. However, in the ensemble 
average, the one body density profiles must be the same and therefore show no signature 
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Figure 4. As in Fig. [3 but now for fluid B with N\ = N 2 = 1000 particles. Initially 
there is no signature of phase separation in the density profiles, but after breaking the 
symmetry in the external potentials, there is. Note the intermediate profiles exhibit two 
fluid-fluid interfaces (rather than one). In the final frame we display a slice (particles 
with — 0.5i?n < z < 0.5i?n) from an equilibrium configuration t < 0. Even though 
the one body density profiles (t — 0) show no signature of phase separation, this 
configuration does. 



of phase separation. Such effects demonstrate an important aspect of the DDFT: this 
is a theory for the ensemble average fluid one body density profiles and therefore, 
as with equilibrium DFT, the fluid density profiles must exhibit the symmetry of the 
confining external potentials. 
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